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ABSTRACT 

Mechanisms regulating the origin of X-rays in YSOs and the correlation with their evolutionary stage are 
under debate. Studies of the X-ray properties in young clusters allow to understand these mechanisms. One 
ideal target for this analysis is the Eagle Nebula (Ml 6), with its central cluster NGC6611. At 1750 pc from 
the Sun, it harbors 93 OB stars, together with a population of low-mass stars from embedded protostars to 
disk-less Class III objects, with age < 3 Myrs. We study an archival 78 ksec Chandra/ ACIS-I observation of 
NGC661 1, and two new 80 ksec observations of the outer region of Ml 6, one centered on the Column V, and 
one on a region of the molecular cloud with ongoing star-formation. We detect 1755 point sources, with 1 183 
candidate cluster members (219 disk-bearing and 964 disk-less). We study the global X-ray properties of M16 
and compare them with those of the Orion Nebula Cluster. We also compare the level of X-ray emission of 
Class II and Class III stars, and analyze the X-ray spectral properties of OB stars. Our study supports the lower 
level of X-ray activity for the disk-bearing stars with respect to the disk-less members. The X-ray Luminosity 
Function (XLF) of M16 is similar to that of Orion, supporting the universality of the XLF in young clusters. 
85% of the O stars of NGC6611 have been detected in X-rays. With only one possible exception, they show 
soft spectra with no hard component, indicating that mechanisms for the production of hard X-ray emission in 
O stars are not operating in NGC 66 1 1 . 

Subject headings: stars: formation, stars: pre-main sequence, X-rays 



1. INTRODUCTION 

Since the discovery of the intense X-ray emis- 
sion from young p re -Main Sequence (PM S) stars 



dFeigelson & DecampliL 11981b iMontmerle et all 11983b . 
X-ray observations of active star-forming regions and young 
galactic clusters become an efficient method to study the star 
formation process and the properties of young stars. In fact, 
the level of X-ray emission in PMS stars, higher than that 
of field main-sequence stars, provides a very efficient means 
of selecting stars associated with star forming regions and 
young clusters. In the last decade a large number of young 
clusters have been observed wi th X-ray telescopes, su ch as 
the Chandra X-ray Telescope dWe isskopf et al., 2002|) and 
XMM-Newton dJansen et all 1200 lb . in order to select their 
young members and study their X-ray activity. To the present 
day, the longest X-ray observations of star-forming regions 
are the 839 k s ec Ch andra Orion Ultradeep Project (COUP, 
iGetman et akl 12005b. the I PS M sec Chandra Cygnus OB2 
Legacy Survey dbrake et afll2009|) and the 1 .6 Msec Ch andra 
Carina Complex Project (CCCP. iTownslev et all 1201 lal) . 

While we have a reasonably detailed comprehension of the 
coronal activity and accretion phenomenon in low-mass pre- 
main sequence stars, and of the mechanisms for X-ray emis- 
sion in massive stars, several key topics are still not com- 
pletely understood. For instance, ho w the presence of a cir - 
cumstellar disk affects X-ray activity dFlaccomio et alir2003l) . 
the nature of the huge spread in X-ray activity obser ved in al- 
most all the young clusters (Feig elson et all 120021) . and the 
importance of the pro posed mechanisms for hard X-ray emis- 
sion in massive stars (Babel & Mo ntmerld 11997b . 

In this paper we address some of these topics by study- 
ing the global X-ray properties of the young stars asso- 
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ciated with the Eagle Nebula (Ml 6) and its central clus- 
ter NGC 6611. Situated in the Sagittarius arm and lying 
in the southern sky, this cluster is 1750 pc away from the 
Sun. It suffers only a low average extinction in the central 
cavity of the molecular cloud clea red by the cluster itself 
(Ay ~ 2.6 m , |Guarcello et al., 2007), but becomes more ab- 
sorbed jn_the_™rrouriding area where the cloud is still dense 
dGuarcello et all I2010bl). Th e central cluster hosts 93 OB 
stars (iHillenbrand et all fT993). including one of the rare can- 
didate 04 stars observed in our g alaxy (W205, with a mass 
of 75 - 80 M , lEvans et alJfeOOSb . the massive binary system 
W175, composed of an Q5V star (lEvans et all 12005) and a 
late-O (likely 08.5) star (Ma rtavan et all 120081) . some Her- 
big Ae/Be sta rs and the two magneti c stars W080 (B1V) a nd 
W601 (B 1 ■5): lMartavan et all d2008l) : lAlecian et all d2008h . 

The region is also rich in Pre-Main Sequence (PMS) stars 
with low- a nd intermediate-mass. In our previou s studies of 
this region dGuarcello et all 120071 120091 I2010allbb . we de- 
tected a total of 1937 candidate PMS stars associated with 
Ml 6. The disk-less stars have been selected from their in- 
tense X-ray emission, and the disk-bearing stars from their 
excesses in infrared bands. The X-ray select ed young stars of 
NGC 6611 have a median age of ~ 1 Myr (Guarcello et al., 
12007b . but in th e whole cloud a star for mation sequence has 
been proposed dGuarcello et all I2010bb . with the oldest star 
formation events to the south-east (~ 3 Myr) and the youngest 
in the north-west (< 1 Myr). Together with these large pop- 
ulations of Class II and Class III o bjects, a significant num - 
ber of embedde d Class 0/1 sources dlndebetouw et all 12007b . 
water masers (Hea lv et all 12004b. and candidate Herbig- 
Haro objects (jMeaburn & Walsh, 1986), have been identi- 
fied in sites of M16 with ongoing star-formation activity, 
such as the structures called "the Pillars of Creation" and 
"Column V" dHester et all [T996t iMcCaughrean & Anderseni 
l2002tlSugitani et all 12007b , the Bright Rimmed Cloud (BRC) 
SFO30 northward of the cluster, and a young cluster in 
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the north-east largely embedded in the mol ecular cloud 
dlndebetouw et all 120071: fGuarcello et ani2009h . 

The Eagle Nebula hosts a large variety of young X-ray emit- 
ting sources whose X-ray properties will be studied in detail 
in this paper. Our scope is to characterize the global X-ray 
properties of the cluster, to study how X-ray activity of clus- 
ter members is affected by the presence of the disks, and fi- 
nally to take advantage of the large massive-star population 
of NGC 661 1 to study X-ray emission in OB stars. The paper 
is organized as follows: in Sect. |2]we review the data analy- 
sis, source detection, the identification of the stellar counter- 
parts of the detected X-ray sources, and the spectral analysis. 
The distributions of plasma temperature and hydrogen column 
density of the Ml 6 stars are analyzed in Sect. [3] The X-ray 
luminosity function of the whole population, as a function of 
the stellar mass and of evolutionary status, is studied in Sect. 
HI in Sect. [5]we study the X-ray emission of massive stars and 
in Sect. [6] we review the X-ray properties of the stellar pop- 
ulation of selected regions of the Eagle Nebula. In Appendix 
|A]we present the X-ray catalog of detected sources. 

2. DATA ANALYSIS AND CATALOG 

2.1. Sources identification and photon extraction 

Fig. Q] shows the image of the Eagle Nebula in the I band 
obtained with the Wiel d Field Imager (W FI) mounted on the 
2.2m telescope at ESO dGuarcello et all [20071) . The contours 
mark the dust emission detected at 8.0 nm with Spitzer/IRAC 
dlndebetouw et all 120071: iGuarcello et all 120091) . marking the 
Pillars of Creation, the Column V and the SFO30 cloud 
(the outer contours delimit the size of the used IRAC im- 
age). The rotated squares represent the fields observed 
with Chandra/ACIS -I. The archival observation of 78 ksec 
dLinskvet all 120071) is centered on NGC 661 1 (the c - field) 
and it includes the Pillars of Creation and the SFO30 cloud; 
the two new observations of 80 ksec (P.I. Guarcello) are cen- 
tered on Column V (the e - field), and on a group of Class I 
and embedded Class II stars (the ne - field) which is younger 
than 1 Myr. 

X-ray image re duction and source dete ction in all the fields 
are described in IGuarcello et all (1201 Obi) , where detected X- 
ray sources (1 158, 363 and 315 in the c— , e—, and ne- fields, 
respectively) with their optical-NIR counterparts have been 
used to classify YSOs in this region. Briefly, for each obser- 
vation a Level 2 event file has been obtained using the CIAO 
4.0 3 tool acis-process-events, retaining only the events incom- 
patible with cosmic rays and removing the 0.5" event position 
randomization added by the standard Chandra data processing 
pipeli ne. Source detection has been performed with PWDe- 
tect 4 dDamiani et all 119971) adopting a threshold correspond- 
ing to 10 spurious detections. The total number of detected 
X-ray sources is 1755 (81 sources are in the overlapping re- 
gion between the c - fie ld and the e - field). As described in 
IGuarcello et al.1 (1201 Obi) , the mass limit of the X-ray sources 
with optical counterparts (for which it was possible to esti- 
mate the mass) is equal to 0.2 M . Fig. Q] shows the photon 
extraction regions of each detected source. The clustering of 
sources in NGC 661 1 is evident, but this image clearly shows 
that also the outer parts of the nebula are well populated by 
young stars. 

Photon extraction has been made using the IDL software 

3 http://cxc/harvard.edu/ciao 

4 http://www.astropa.unipa.it/progetti_ricerca/PWDetect 



ACIS Extract 5 (AE, iBroos et all 120101) . which uses TARA 6 , 
CIAO, FTOOLS 7 and MARX 8 packages. AE calculates the 
Point Spread Function (PSF) for each source, using it to de- 
fine a photon extraction region as the region encompassing 
90% of the PSF evaluated at \ A9keV. Background events 
have been extracted for each source in circular annuli centered 
on the sources, applying a mask for source counts defined in 
two iterations (first a circular area covering 99% of the local 
PSF and then with a more accurate mask region). The extrac- 
tion regions of crowded sources are reduced in order to not 
overlap each other, down to 40% of the local PSF in the most 
crowded regions. In this first phase, AE uses the source posi- 
tions provided by the source detection algorithm (PWDetect 
in our case). Once both the PSF and the background for each 
source are evaluated, AE computes new source positions by 
correlating the source images with the local PSF model. The 
user has the possibility to use these new positions or retain the 
set of positions provided by the source detections algorithm. 
As suggested by the AE manual, we updated the positions of 
the sources observed with an off-axis larger than 5' (a fur- 
ther coordinates correction of A a = -0.18" and = -0.05" 
has been applied to cross-correlate the X-ray catalog with the 
optical-infrared catalog, see Sect. |2.2| i. AE, then, repeats the 
computation of the PSF and the background and the photon 
extraction using the new positions. For each source, AE pro- 
vides information such as the net observed counts, the light 
curves and a probability that the observed light curve is con- 
stant based on a Kormogorov-Smimov test; the most relevant 
are summarized in the electronic catalog described in Ap- 
pendixlAl 

2.2. Stellar counterparts of the X-ray sources 

Detec ted X-ray sources are classified in IGuarcello et al.1 
(2010b) according to the properties of their optical/infrared 
counterparts. Among the 834 sources with disks, selected by 
their NIR excesses with respect to the expected photospheric 
emission, 219 have an X-ray counterpart and are classified as 
"disk-bearing members" . We classified as "disk-less mem- 
ber" 964 X-ray sources that are younger than 10 Myr and 
with an extinction higher than 2.6™, according to their po- 
sition in the optical color-magnitude diagram. A total of 76 
sources are older than 10 Myr and/or have extinction smaller 
than 2.6'". These sources are most likely foreground main 
sequence stars detected in X-rays, and then they have been 
classified as "foreground sources". A total of 504 X-ray 
so urces are not id entified with a known stellar counterpart. 
In IGuarcello et al.1 d2010bl) we compared our count-rate limit 
(6.25 X 10~ 5 counts js), properly converted into a limit flux 
in the 0.5 - lOfceV energy band, with the LogN vs. Log S 
di stribution of extrag alactic sources in the ELAIS field shown 
in iPuccetti et al.l d2006l) . estimating an upper limit of the ob- 
served extragalactic X-ray sources equal to 193.5. We then 
expect that about 300 X-ray sources without optical-infrared 
counterparts can be very embedded sources associated with 
M16 (a significant fraction of them lie in the regions with high 
extinction, such as the ne- field) or background sources lying 
in our Galaxy. Since we cannot discern among these hypothe- 
ses, we simply classified them as "background sources " and 
we have not considered them in our study of the global X-ray 

5 http://www2.astro.psu.edu/xray/docs/TARA/ae_users_guide.html 

6 http://www.astro.psu.edu/xray/docs/TARA 

7 http://heasarc.gsfc.nasa.gov/docs/software/ftools 

8 http://space.mit.edu/CXC/MARX/ 
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Fig. 1. — Image in I band of M16, with the contours encompassing the emission at 8.0jum obtained from Spitzer/IRAC. The rotated squares are the 
Chandra/ACIS-I fields of the observations analyzed in this work. The black circles mark the positions of detected X-ray sources. The radius of each circle 
is proportional to the PSF area at sources position. The positions of Pillars of Creation, Column V, the SFO30 cloud, the North-East embedded cluster, and 
NGC 661 1 are also shown. 



properties of the Ml 6 population. 

Fig. [2] shows the spatial distributions of these four classes 
of X-ray sources. The disk-bearing stars are clearly clumped 
in NGC 6611, toward the Column V and in the embedded 
north-east cluster. The disk-less members are more numerous 
and distributed across the whole cloud. Also the background 
sources are sparse in the field, but the high concentration of 
them in the regions with the highest extinction and embedded 
ongoing star formation activity (the ne - field, eastward the 
Column V, and in the tip of the Pillars of Creations and Col- 
umn V) is consistent with the idea that most of these sources 
are very embedded young stars associated with M16. 

2.3. Spectral analysis: spectral fitting 

The determination of the absorption corrected X-ray lu- 
minosity (Lx), as well as the plasma temperature (kT) and 
hydrogen column density (Nh), requires the analysis of the 
X-ray spectra. AE provides both the source and back- 



ground spectra, the redistribution matrix files (RMF) and the 
ancillary response files (ARF). We fit the observed spec- 
tra with thermal plasma (with both one and two tempera- 
tures) and power-law models. We u se the APEC ioniz ation- 
equilibrium thermal plasma code (Smith et al., 200 fl), as- 
sumin g the sub-solar elemental abundances of Maggio et al. 
(2007) ; the absorption was trea ted using the WABS model 
(Morri son & McCammonl Il983l) . The one-temperature ther- 
mal model was applied to all the sources with more than 25 
counts; while the two temperature thermal model was applied 
to each source with more than 80 counts. The power-law 
model has been applied to those sources with hard spectra for 
which the best-fit thermal model predicts a plasma tempera- 
ture kT > 5 keV. In order to avoid false convergences due 
to local minima in the x 2 space, we applied different sets of 
initial values to each model. When more than one model has 
been used for a given source, we choose the best model by the 
X 1 probability and visual inspection of the spectrum. 
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In the c - field we were able to obtain a good spectral fit 
for 441 sources, 1 14 in the e - field and 60 in the ne - field 
(out of a total of 1 158, 363 and 315 sources, respectively), for 
a total of 575 sources with a good spectral fit (40 sources in 
the overlapping regions). For 33 sources the best fit model 
is a two temperature thermal plasma model. For a total of 
20 sources the best fit model is a power-law spectrum: 15 of 
them are classified as "background sources ", the remaining 5 
as "disk-less members". 

2.4. Spectral analysis: energy quantiles 

Spectral analysis cannot be performed for faint sources. For 
them, an estimate of the X-ray luminosity can be provided 
by th e analysis of the photon energy quantiles (Hon g et ail 
2004). This method is similar to the analysis of the X-ray 
hardness ratio, adopted in several works, but it does not suffer 
the drawback of the choice of the energy bands us ed to define 
the "X -ray colors". Adopting the definition by Hong et al. 
(2004), if x% is the percentage of the total counts with energy 
below the value E x %, the quantile Q x is defined as: 

Qx = ^Z^ (1) 

E up lo 

where E/ and E up are the lower and upper limits of the en- 
ergy range, respectively 0.5 keV and 8.0 keV in our case. We 
used the median energy (i.e. the 50% quartile) and the 25% 



and 75% quartiles, which are combined in the following in- 
dependent variables: log(Q 50 / (1 - g 50 )) and 3 x {Q 2 slQi5)- 
In order to estimate the values of Nh and kT for each source, 
the observed quantile variables have been compared (by in- 
terpolation) with a theoretical grid of thermal models with 
sub-solar abundances and a gri d of power-law mo dels with 
spectral indices typical of AGN (Brandt et all 120011) . 

Fig. [3]shows the thermal and power-law grids and the quan- 
tiles of the X-ray sources (only the sources with a significance 
larger than 2 are shown; see Appendix[A]) for the four groups 
of sources defined in Sect. 12.21 The thermal models have a 
plasma temperature ranging from 0.3 keV to lOkeV and hy- 
drogen column density from 0.01 x 10 22 cm~ 2 to 6 x 10 22 cm" 2 ; 
the power laws have the same hydrogen column density range 
and a photon index ranging from 0.1 to 2.0. The actual grids 
used for the interpolations are finer than those shown in Fig. 

El 

The disk-bearing and disk-less sources (upper panels) are 
mostly concentrated in the region with 1 .2 keV < kT < 4 keV 
and from 0.1 x 10 22 cm' 2 < N H < 1.1 x 10 22 cm 2 , cor- 
responding to a wide range of visual extinction. Both the 
background and the foreground sources are mostly outside 
the thermal grid. The distribution of the latter suggest val- 
ues of Nh smaller than those used for the grid. The former 
mostly show harder spectra and, in most cases compatible 
with power-law models. However, from the grid it is not pos- 
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Fig. 3. — Quantiles of the four classes of X-ray sources, with overplotted the grids derived from thermal one-temperature spectra (blue lines) and power-law 
spectra (red lines). Some of the values of Nh and kT for the thermal models are shown. The power law models have indices ranging from to 2.0 and the same 
Nfj of the thermal models. 



sible to discriminate the real extragalactic sources from the 
embedded cluster members, given the degeneracy between 
very hard thermal spectra and power laws. 

A close inspection of Fig. [3] shows the presence of 8 back- 
ground sources in the area of the soft thermal spectra cor- 
responding with kT < 0.5 keV. Even if they are not very 
bright (the median net counts of this sample is 19), the prob- 
ability associated with the hypothesis that they are spurious 
detections is negligible (see the PBS parameter, explained in 
Appendix [A}. Five of these sources lie in the c - field and 
have a median photon energy between 1.3 keV and 1.8 keV. 
It is not obvious what these sources are. One intriguing pos- 
sibility is that some of them are compact objects associated 
with Ml 6, which are remnants of the explosion of a massive 
star. The possibility that Ml 6 h osted a supernova e xplosion in 
the past has been suggested by Flagev e t alj (1201 lb to explain 
the hot dust (~ 70^) they observed in the cloud. However, 
NGC 6611 is not old enough to firmly claim that the most 
massive members already exploded as supemovae, and the 
median photons energy of these 8 sources is higher than that 
of the candi date neutron stars recen tly identified in the Carina 
complex bv lTownslev et al.l (1201 lbl) . 



In total, 846 sources lie inside the thermal grid. For these 
sources the individual values of Nh and kT have been calcu- 
lated by interpolating their quantile variables in the thermal 
grid, and then their X-ray luminosi ty, Lx, using the di stance 
adopted for the cluster (1750 pc, Guar cello et al.ll2007h and 
correcting for the individual absorption. To estimate the un- 
certainties in the X-ray stellar parameters from the grid, we 
calculated the errors in their photons ene rgy quantiles, fol- 
lowing the procedure of Hong et al. (2004) and propagated to 
the independent quantile quantities. Then, we repeated the 
interpolation with the grid four times, each time with a dif- 
ferent combination of sum and subtraction of the errors to the 
quantile quantities. In this way, we had, for each star and for 
each parameter, 4 slightly different values. For each nomi- 
nal value (i.e. those derived from the grid without adding or 
subtracting the errors) of the three quantities, the upper er- 
ror has been defined as the difference between the maximum 
of this 4 measures and the nominal value, the lower error as 
the difference with the minimum value. The median error for 
the three parameters are <tn h = 0.7 cm" 2 , <Jk T = QAkeV, and 
o~h g (L x ) = 0.3 erg/s. 

A total of 312 sources lie inside to the power-law grid, and, 
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as shown in Fig. [3] 82% of them are background sources 
(with a wide range of hydrogen column density) and disk-less 
stars. There are 346 sources lying outside both grids. Most 
of them have 3Q?s /Q7S smaller than those of the grid models. 
iHong et al] (120041) suggested that low significance sources can 
fall in the area of the diagram with 3 625/675 < 0.8. In fact 
the sources in this region of the diagram have a median of 12 
counts, indicating they mostly fall in this region because of 
low-counts. 

2.5. Upper limits 

The main aim of this paper includes the comparison be- 
tween the X-ray properties of disk-less and disk-bearing stars, 
the study of the global X-ray properties of the M16 popula- 
tion, and their comparison with those of the Orion members. 
A reliable comparison requires the knowledge of the upper 
limit of X-ray luminosity for the undetected members. The 
problem of the incompleteness of the sample of Class III ob- 
jects of M16 will be discussed later, since we do not know the 
X-ray undetected Class III stars. This is not the case for the 
disk-bearing members, since we have detected in X-rays 219 
sources with infrared excesses out of a total of 834 previously 
identified. The upper limit of the photon count rate at the po- 
sitions of the undetected disk-bearing stars is calculated with 
PWDetect, adopting the same threshold significance used for 
source detection. We have converted these upper limit count 
rates to X-ray luminosity using as conversion factor the me- 
dian of the ratio between the X-ray luminosities and photon 
count rates for the detected members. 

3. PLASMA TEMPERATURE AND HYDROGEN 
COLUMN DENSITY OF M16 MEMBERS 

In the following sections we will address several facets of 
the X-ray properties of M16 members, starting from the distri- 
butions of plasma temperatures and absorption. As explained 
in the previous sections, values of Nh and plasma temperature 
are associated to the X-ray source from the spectral fitting 
or quantile analysis. When a good fit with a thermal model 
is available, these values and their uncertainties are adopted 
from the best fit model. In the other cases, they are taken from 
the interpolation in the thermal grid (shown in Fig. |3) for the 
sources which lie inside the grid. For those sources without a 
good spectral fit with a thermal model (IT or 2T) and whose 
position in the diagrams in Fig. [3]is outside the thermal grid, 
no values of Nh and plasma temperature are provided. 

Fig. H] shows the distributions of Nh, kT and median pho- 
ton energy of the X-ray sources in M16 (both disk-less and 
disk-bearing sources). The Nh distribution covers a large 
range of values, with a median value corresponding to Ay = 
3.6 m , which is larger th an the average extinction estimated by 
iGuarcello et al.l (120071) from the analysis of the optical dia- 
grams (Ay = 2.7™). We will show in Sect. [6]that the discrep- 
ancy between these two estimates is reduced by considering 
the stars in the central cavity cleared by NGC 6611. Part of 
the large spread of the Nh distribution is purely statistical, 
but since about 55% of sources with lowest Nh are clustered 
in the cavity cleared by NGC 661 1 or southw ard (where the 
extinction is lower, see IGuarcello et al.ll2010bl) . while a sim- 
ilar fraction of the sources with highest Nh mostly lie in the 
trunks, in the SFO30 cloud and in the north-east, part of the 
spread in the Nh distribution is due to the strong differential 
reddening across the nebula. 

The median values of both kT and median photon energy 
are similar to those observed in other star-forming regions (for 
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instance in Orio n, Feigel son etail 120021 and in NGC 1893, 
iCaramazza et al.L 1201 lh . Thetwo populations of sources with 
kT < \.\AkeV and kT > 4.09 keV (which are the quantiles 
at 80% and 20% of the kT distribution), have different frac- 
tions of variable sources, suggesting that part of the spread 
observed in the kT distribution can be due to X-ray variability, 
with some sources observed during periods of quiescent emis- 
sion and others during flares. In fact, by adopting a threshold 
for variability of P KS < 0.005 (see Appendix lAl. 42% of the 
sources in the high-temperature tail of the kT distribution can 
be variable, while among those with kT < 1.14 keV only 14% 
can be variable. This suggests a larger contribution from flares 
in the X-ray emission of the sources with highest plasma tem- 
perature. Besides, the presence of sources with high plasma 
temperature [kT > 4.09) among those with constant light 
curve and low luminosity supports the hypothesis that "quies- 
cent" X-ray emission from magnetically a ctive stars can arise 
from the merging of co ronal microflares (iDrake et al.l |2000[ 
ICaramazza et alll2007h . 

Fig. [5] shows the distributions of Nh, plasma temperature 
and median photon energy separately for sources in the c - 
field (left panel), e— field (central panel), and ne- field (right 
panel). The distributions for the populations of the c - field 
and e - field are similar, reflecting t he similar age and extinc - 
tion in these two fields as found by IGuarcello et al.l (|2010b). 
The largest differences arise with the ne - field, where the 
median hydrogen column density is Nh = 1.4 x \0 22 cmT 2 , 
corresponding to Ay = 7.98"', which is higher than the extinc- 
tion and median Nh found in the other two fields (3.2 m ). This 
mostly explains the lack of sources with soft spectra in the 
ne - field compared with the other two fields, and the higher 
kT median value (2.86 keV versus 2.01 keV and 2.12 keV). 

4. X-RAY LUMINOSITY OF M16 MEMBERS 
4. 1 . The X-ray Luminosity Function 
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Fig. 5. — Distributions of the Nh, kT and median energy of the X-ray sources in M16 falling in the three observed fields. The dotted vertical lines mark the 
25%, 50%, and 75% quartiles of each distribution. 



X-ray luminosities and their uncertainty have been obtained 
from the spectral or quantile analysis with the same criteria 
adopted for plasma temperature and hydrogen column (Sect. 
|3). In addition, for the sources without a good spectral fit and 
that are outside the thermal grid in the diagrams in Fig. [3] 
they have been inferred directly from the count rates, using as 
conversion factor the median Lx/ '{count .rate) ratio found for 
the other sources (which ranges from 7 x 10 33 to 1 x 10 34 in 
the three fields). 

The total luminosity of Ml 6 is equal to log(Lx) = 
34.90 erg ■ s~ l , with an individual source median value of 
log(Lx) = 30 A^ A 6 erg ■ s~ l . In this estimate we excluded 
a candidate disk-less member with very high luminosity 
log(Lx) - 34.1^ erg ■ s -1 . This luminosity, in fact, is 
apparently too high to be correct, since it is two orders of 
magnitude higher than the Lx of the second most luminous 
star in our sample. Its spectrum has been fitted with a 2 
temperature thermal model, with a soft primary component 
(kT[ = 0.11 + 0.01 keV) and a harder secondary component 
with kTi = 1.7 + 0.6 keV, but the emission measure of the 
hard component is 10 5 times fainter than that of the soft com- 
ponent. One possibility is that this is a foreground source with 
an overestimated distance, or a background High Mass X-ray 
Binary, but both the N H value (equal to 4.45 x 10 22 cm' 1 ) and 
its position in the optical color-magnitude diagrams suggest a 



young and embedded star lying in the BRC-SFO30. 

Fig. [6] compares the X-ray luminosity calculated from 
the quantile analysis with that from the spectral fits, for the 
sources in our sample with a good significance and a spectral 
fit with good^- 2 . The estimate of Lx from the quantile anal- 
ysis agrees with that from the spectral analysis for most of 
the stars in Fig. [6] with the exception of few outliers. These 
outliers mostly lie in the low kT-high Nh region of the grid, 
where the low sensitivity of the ACIS detector to the soft part 
of the spectrum result in large uncertainty in estimates of Nh 
and kT. 

One of the aims of this paper is to verify whether the X-ray 
Luminosity Function (XLF) of M16 confirms the universality 
of the XLFs in young clusters younger than few Myrs, as a 
consequence of the very slow decline of the X-ray luminosity 
for PMS stars younger than 10 Myrs (Preibisch & Feigelson, 
120051) . Several authors compared the XLFs of young clusters 
with that observed in the Orion Nebula Cluster (ONC) in the 
COUP survey, which is the most complete X-r ay observation 
of a y oung cluster, complete down to 0.1 M Q (Getman et al., 
2005). A good agreement was f ound i n IC348 and NGC 133 3 
(2-3 Myrs, iFeigelson etafl 120051: iWinston et all 120101) , 
M17 (~ 1 Myr, iBroos et al.L 12007b. NGC 6357 (~ 2 Myrs 
Wang et all |2007|), NGC 2244 (2 - 3 Myrs, IWang et all 
2008), and Cyg OB2 (3-5 Myrs, Wright et al., |20Jj0j). In 
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Fig. 6. — Values of the X-ray luminosity (in erg ■ s ) calculated from the 
quantile analysis vs. the values obtained from spectral fitting for the X-ray 
sources in our sample. 




Fig. 7. — X-ray Luminosity Functions for M16 and Orion members. The 
dashed histogram is the Orion XLF scaled up to the Ml 6 XLF. The con- 
tinuous and dotted lines shows the slope of the linear fits of the XLFs for 
log(Lx) > 30 erg ■ s' 1 . 



this paper we adopt the same approach, comparing the XLF 
of Ml 6, where the X-ray emitting popula t ion ha s an age rang- 
ing from < 1 - 3 Myrs dGuarcello et all 120071) . with that of 
the ONC. 

Fig. |7]shows the XLF we obtained for M16, together with 
that of the COUP sources in Orion. The latter is almost con- 
stant for log (Lx) < 30 erg ■ while the M16 XLF drops be- 
low this limit, which turns out to be the completeness limit of 
our observations. For luminosities higher than the complete- 
ness limit, the XLF of Ml 6 can be fitted with a power law 
with index F = -0.85 + 0.09, whic h is consistent w i th that of 
the Orion XLF (r = -0.93 + 0.08. iFeieelson et all 120051) in 
the same luminosity range. In this calculation, we have not 
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Fig. 8. — log (Lx) vs mass for the detected Class III members in M16 (dia- 
monds) and Orion (dots). 



considered the OB stars detected in X-rays, since the Orion 
OB population is significantly different from that of M16. 

The comparison between the XLFs of Orion and M16 can 
help us to estimate the level of incompleteness of our sam- 
ple of members down to the subsolar limit. To match the 
Orion XLF with the M16 XLF, we have to scale the former 
by 1.7. The tot al Orion population with M > 0.1 M Q counts 
1600 members (Get man et alll2005l) . The match between the 
two XLFs suggest a total M16 population of 2700 members, 
larger than the number of members we identified in our previ- 
ous work (1907 in total, incomplete a low mass regime). This 
is a rough estimate of the still non identified M16 population, 
since it depends on several properties of M16 and ONC, such 
as the physical extent and absorption characteristics in the re- 
spective sky areas, that are not compared properly here. 

4.2. X-ray luminosity and stellar mass 

In young PMS stars, Lx is generally observed to increase 
with the stellar mass, due to the fact that the dynamo mecha- 
nism in PMS stars is in the saturated regime, which means 
that it is independent of the rotational velocity and it de- 
pends only on the ste l lar m ass and bolometric luminosity 
dPreibisch & Feigelsonll200"5l) . To study the Lx vs. mass rela- 
tion in Ml 6, we estimated the masses of M16 members with 
optical and X-ray detections, first obtaining their intrinsic op- 
ti cal photometry using th e extinction map of M16 obtained 
in lGuarcello et al.l (120 10b). then finding both stellar mass and 
age by interpolating the dereddened V and V - 1 with t he col - 
ors expected from the grid of isochrones of lSiess et al.l (12000). 
In this study, we focused only on the Class III stars, since the 
mass estimate from the optical photometry in Class II objects 
can be affected by the presence of the disk ( Gu arcello et all 
120103) . 

Fig. [8] shows the Lx vs. mass sca t ter pl ot for Class III 
stars in M16 and Orion (Getman et al., 2005). In each mass 
bin, there is a large spread of Lx of about 1.5-2 orders of 
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magnitude. This spread is mostly statis tical, but also binarity 
can play a role (Feigelso n et al.L 1200 21). The distribution of 
the low-mass stars in M16 (i.e. with M < 0.8 M Q ) follows the 
Lx vs. mass relation valid in Orion. In fact, a linear fit in this 
range of mass gives: 

log(Lx) = (30.9 ± 0.1) + (1.1 + 0.3) ■ log(M/M B ) (2) 

with a slope consistent with that found in Orion 
( Preibi sch & Feigelsonl 120051) . Given the consistency of the 
two XLFs and assuming the universality of the IMF, the fact 
that the low-mass population of Ml 6 and ONC follow the 
same Lx vs. mass relation is not surprising. 

In Fig. [8]the distribution of stars more massive than 0.8 M 
is totally different, being flatter (with a slope of 0.4 + 0.2) 
than in the low-mass regime. Similar behavior has been ob- 
served in other young clusters. Our data do not allow us to 
understand the nature of the flattening of the Lx vs. mass dis- 
tribution observed at about 1 M . The most li kely explanation 
is that it is due to changes in st ellar interior. Feigelso n et al.1 
(120031) . following the models of lPalla & Stahleri (119931) . have 
shown as in the mass range between 2.4 < M/M Q < 3.9, 
very young PMS stars already developed a radiative core, sur- 
rounded by a convective mantle, where deuterium burning 
takes place. The efficiency of the magnetic dynamo, and the 
X-ray activity level of the star, strongly depends on the pre- 
cise boundary between these two regions, and this boundary 
changes with the age and accretion history of the PMS phase. 
For even more massive stars, the mechanism for the produc- 
tion of X-rays changes drastically. 

4.3. X-ray activity in Class III and Class II members 

It is still debated whether the X-ray activity in PMS 
stars can be affected by the presence of a circumstellar 
disk. Several authors have compared the X-ray luminos- 
ity of Class III and Class II members of young clusters, 

and the general pictur e remains uncle ar. In some cases, 

such as Chama l eon I dFeigelson et al.L 119931). o Ophiuchi 
(Casa nova etafl 11995b and IC348 ([Preibisch & Zinneckerl 
120021) . the same level of X-ray activity in disk-less and 
disk-bearing members has been found; while in other clus- 
ters, s uch as the Ta urus- Auriga region (IStelzer & Neuhauserl 
l200l. NGC 1893 (iCaramazza et all 120121). and NGC 2264 
(Flaccomio et al. , 2006), the X-ray luminosity of Class II ob- 
jects is significantly lower than that in Class III stars. In 
the ONC the first attempt to find any difference in X-ray 
act ivity betwee n these two classes of PMS stars was made 
bv lGagne et al.l (119951) , who found no significant difference. 
Later. IFlaccomi oet al.l (12003 ) discovered a different emission 
level in accreting and non accreting stars in Orion. This re- 
sult was reinforced by Preibisch & Fei gelsonl (12005b . taking 
advantage of the completeness of the member selection with 
the deep COUP data. These authors found no difference in 
X-ray emission between stars with and without infrared ex- 
cesses, detected with a diagnostic based on the L photometric 
band, but they found a significant difference between non ac- 
creting and accreting stars, selected using the Call line emis- 
sion. Their finding suggests that the effects on the X-ray 
activity in disk-bearing stars are not due to the presence of 
the disk itself, but to active ac cretion. In fact, in the ONC 
study of P risinzano et al.l d2008l) . where the disk-bearing stars 
have been selected by means of their infrared excesses, the 
level of X-ray emission in Class II and Class III stars is only 
marginally different. 



The present study of the stellar population of M16 can help 
to shed some light on the effects of the circumstellar disk in 
the X-ray emission in young stars, taking advantage of the ac- 
curate selection of both Class II and Class III members. How- 
ever, a reliable comparison between the X-ray luminosity dis- 
tributions of Class II and Class III objects can be attempted 
only if the two samples are complete. Class II sample (219 
stars) can be easily completed using the u pper limits of Lx of 
the undetected Class II objects (see Sect. 12.5b . These upper 
limits, together with the X-ray luminosity of detected Class II 
and Class III objects, are shown in the left panel of Fig. [9] as 
a function of stellar mass. 

We cannot complete the Class III sample with the same ap- 
proach adopted for Class II stars, since we have not a rep- 
resentative sample defined independently from X-ray infor- 
mation. However, we can estimate the fraction of Class III 
members with X-ray emission below our completeness limit 
using the COUP data of Orion. Our sample of Class III ob- 
j ects, as well as the Cla ss II, is more massive than 0.2 M 
dGuarcello et al.L 1201 Obi) , and in Orion 15% of the Class III 
objects with M > 0.2 M have log(L x ) < 30 erg ■ s -1 . Given 
the similarities between the X-ray properties of the two clus- 
ters, we can suppose that this is valid also for Ml 6, i.e. that 
the undetected Class III stars more massive than 0.2 M e are 
15% of the detections (124 sources). Regarding their upper 
Lx limit, since our aim is to verify whether the X-ray activ- 
ity in Class III sources is higher than that in Class II objects, 
we can consider the most conservative case, resulting in the 
lowest Lx distribution for Class III stars: that the upper limit 
of the X-ray luminosity of the missing Class III stars is equal 
to the l owest Lx measured f or Class III objects in the COUP 
survey (Getman et al., 2005), equal to log (Lx) = 27.89 erg/s. 
We then calculated the XLFs for the "completed" samples of 
Class III (detected Class III sources plus 124 sources with 
log (Lx) = 27.89 erg/s) and Class II objects (detected Class II 
objects plus the upper limits of the undetected Class II ob- 
jects.), which are shown in in the right panel of Fig. [9] The 
emission level of the completed sample of Class II objects is 
evidently lower than that of the Class III, and we h ave ver- 
ified with the ASURV 9 dFeigelson & NelsonL 119851) statisti- 
cal package that the probability that the two XLFs are drawn 
from the same parent distribution is null. In the left panel, 
a few sources with log(Lx) = 29.5 erg/s are present (fainter 
than our detection limit). They are low significance sources, 
whose luminosity has been inferred directly from their count- 
rates. In conclusion, in M16 the level of X-ray emission is 
higher in Class III objects than in disk-bearing stars, support- 
ing the idea that the presence of a disk results in a lower X-ray 
activity. We cannot discern between accreting and non accret- 
ing disks, but it can be expected th at, since the average a ge 
of the M 16 population is ~ 1 Myrs dGuarcello et al.L 120071) . a 
large fraction of stars with disks is still actively accreting. 

There is more than one hypothesis to explain the observed 
lower X-ray emission level in Class II than Class III stars. 
The accreting material, for instance, could be responsible for 
a higher absorption in disk-bearing objects. In fact, the cor- 
rection we made for the extinction can be not completely ap- 
propriate for disk-bearing stars. For these objects, it would be 
more appropriate the use of an extinction model which takes 
into account the presence of both interstellar and circumstel- 
lar obscuring material (Flaccomio et al., in preparation). It 
is also possible that the accreting material deforms the large 

9 Astronomy Survival Analysis, http://www.astro.psu.edu/statcodes 
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Fig. 9. — Left panel: log(Lx) vs. M/M Q for the Class III (crosses), Class II (filled dots) and upper limits (arrows) in M16. Right panel: X-ray luminosity 
function for completed Class III and Class II sources with mass in the 0.2 < M/M Q < 7 range. 



scale stellar magnetic field ( Roma no va et al.L 120041) or that it 
increases the gas density in the magnetic loops around the 
stars resulting in a less efficien t heating of the materia l by th e 
energy released during flares (Preibisch & Feigelson, 2005), 
even if this would require a fraction of the stellar surface cov- 
ered by hot spots larger than the typical values observed in 
CTTS dMuzerolle et al.L 120011) . The last possibility is that 
the accretion phenome non changes the internal struct ure of 
the star, as pointed by IWuchterl & Tscharnuterl (12003b . who 
found that even with small accretion rates the central star is 
no longer fully convective. A different hypothesis is that the 
presence of the accreting inner disk is not the cause of the 
lo wer X-ray ac t ivity, but a consequence, since, as discovered 
by iDrakeet al. (2009), a high X-ray flux emitted by the cen- 
tral star can enhance the photoevaporation of the disk. 

Only one Class I object has been detected in X-rays. Its po- 
sition is a = 20 : 33 : 51.237 and 6 = +41 : 25 : 10.68, and 
it has only IRAC and Chandra detections with a high X-ray 
luminosity log(Lx) = 31. 4-9 erg ■ s _1 derived from the spec- 
tral fit with a 1 temperature thermal model. This source lies 
in the embedded cluster in the ne - field, and it is not surpris- 
ing that it has a very high extinction (Nh — 48.0 x 10 22 cmr 2 , 
corresponding to Ay = 60. 8 m ). The X-ray luminosity of this 
source is about half order of magni tude higher than that o f 
the brightest Class I objects in Orion (Prisinz ano et al.L 2008), 
but a simila r embe dded object has been identified in Ml 7 
dBroos et al.L 12007b . However, neither the IRAC or the X- 
ray data allow us to discard the possibility that this source is a 
background AGN, whose locus in the IRAC color-color dia- 
gram (used to classify the selected disk-bearing stars) roughly 
coincides with that of Class I objects. The fit of the source 
spectrum with a power law model is statistically acceptable 
and the photon index we obtained is equal to 2.03, compatible 
with typical values of AGN spectra. The only valid argument 
against the extragalactic nature of this source is its spatial cor- 
respondence with the embedded NE cluster. 

5. X-RAYS EMISSION IN HIGH MASS STARS 



NGC 6611 is a young massive cluster comprising several 
massive stars. In the ACIS FoVs fall 13 O stars and 80 B 
stars, with one of the few 04-05 stars known in our Galaxy 
with a mass of 75 - 80 M Q (W205), and two known O binary 
systems (W175, an 08.5V+Q5 V system. jBosch et all [19991 
and W205 itself, Q4V+Q7.5V. ISana et all 12009b . In our ob- 
servations, we detected 85% of the O stars (11/13), 48% of 
the B0-B2 stars (24/54), and 19% of the B3-B9 stars (7/36), 
among which 2 giant stars, accounting for a good sample of 
stars to study the mechanisms for the production of X-rays in 
young mas sive stars. In th i s secti on we use the stellar nomen- 
clature by lde Winter et all (119971) . 

The X-ray properties of the detected OB stars in NGC 661 1 
are summarized in Table Q] For those stars observed in two 
observations (in all the cases in the c - field and e - field), 
the table shows two values of Nh, kT, and Lx, corresponding 
to the c - field and e - field. The variability diagnosis comes 
from the KS test on the light curve provided by AE (see Ap- 
pendix | A| for deta il s); the s pectral class i fication is taken from 
iHillenbrand et all d 19931): Evans et all d2005l) : iDufton et all 
( 2006): iMartavan et al.l (120081). w h ile the binarity i n forma - 
tion is from iDuchene et al] (I2001h: iHerbig & Dahml (120011) : 
iGvaramadze & Bomansl d2008l) : ISana et al] (120091) . 

5.1. X-ray emission from the O stars 

Generally, O and late B stars are sources of soft 
(kT ~ 0.5 keV) X-rays, thought to be produced in a myr- 
iad of small shocks in their radiatively accelerated winds 
dOwocki & Cohe n. 1999). This emission is usually slowly 
variable, and, since the wind intensity scales with L/ m i, X- 
ray emission also scales with Lb \, usu ally approximately as 
L x = 1Q- 7 xL/w = 10 3 1 - 10 33 erg • s 1 (iHarnden et all [1971 
Pallav iciiii et all 119811) . However, a large scatter is usually 
observed, mostly due to a wide range of shock filling factors 
and binarity. Usually, a thermal emission spectrum observed 
in an OB star with kT = 0.5 - 0.7 keV and with constant or 
slightly variable light curve is a signature that the X-ray emis- 
sion is produced by the self-shocked wind. In some cases, 
however, also an intermediate/hard component in the X-ray 



TABLE 1 

X-RAY PROPERTIES OF THE OB STARS. 



ID 


net counts 


results* 


N H (xW 22 cm- 2 )" 


kT (keV)" 


kT" 


log (Lx) {erg ■ s ')** 


log (Lboll L sun ) 


Variable*** 


spectral type 


binar 


W205 


754.7 


F 


0.54 


0.63 




32.17 


5.58 


Y 


04III,04V,05 


Y 


W175 


1258.0 


F 


0.78 


0.37 


3.62 


32.65 


5.49 


Y 


05V+08.5V 


Y 


W197 


1041.0 


O 


0.05 


0.93 




31.47 


5.23 


Y 


07V 


Y 


W222 


176.2 


F 


1.69 


0.36 




32.09 


5.14 


N 


07V07III 




W246 


272.2 


F 


0.67 


0.53 




31.70 


5.62 


N 


0711 




W401 


413.3 


F 


0.15 


0.54 




31.26 


4.86 


N 


08.5V 




W161 


208.6 


o 


1.24 


0.34 




32.33 


4.86 


N 


08.5V 


? 


W166 


151.4 


F 


0.45 


0.54 




31.06 


4.82 


N 


08.5,09V 




W314 


99.6 


F 


0.19 


0.62 




30.64 


4.52 


N 


O9V,B0V 


Y 


W280 


60.9 


F 


0.51 


0.59 




31.05 


4.68 


Y 


09.5Vn 




W367 


119.6 


F 


0.01 


0.44 




30.59 


5.15 


N 


09.7III,09.5V 




W188 


36.6 


F 


0.74 


0.97 




30.55 


4.26 


N 


B0V 


Y 


W259 


21.5 


N 


0.14 


10.0 




30.06 


4.04 


Y 


B0.5V 




W150 


19.91 


Q 


0.05 


1.01 




29.62 


4.04 


Y 


B0.5V 




W469 


12.6 


Q 


1.12 


0.70 




30.51 


4.04 


? 


B0.5Vn 


? 


W351 


15.7 


Q 


0.05 


2.90 




29.76 


3.82 


N 


B1V 




W343 


33.6, 30.5 


F 


0.09, 0.67 


22.3, 2.62 




30.24, 30.30 


3.82 


Y 


B1V 


Y 


W25 


65.5 


F 


0.83 


1.49 




30.82 


3.93 


N 


B1V,B().5V 


Y 


W125 


29.6 


Q 


0.33 


1.67 




30.17 


3.82 


Y 


B1.5V,B1V 


Y 


W207 


196.0 


F 


0.51 


1.95 




31.07 


3.82 


Y 


B1V 




W231 


156.3 


F 


0.49 


2.97 




30.99 


3.82 


Y 


B1V 




W254 


3.58 


Q 


0.19 


0.30 




29.09 


3.82 


? 


B1V 


Y 


W444 


17.7 


Q 


0.09 


2.08 




30.02 


3.63 


? 


B1V, B1.5V 


? 


W296 


42.3 


F 


0.17 


1.88 




30.23 


3.57 


N 


B1.5V 




W601 


31.3 


F 


0.11 


1.72 




30.22 


3.57 


Y 


B1.5V 




W421 


93.5, 118.6 


F 


0.48, 1.67 


2.64, 1.15 




30.81, 31.30 


3.57 


N 


B1.5:V 




W80 


42.1 


F 


0.09 


4.14 




30.33 


3.57 


N 


B1VB2V 




W300 


14.5 


Q 


0.05 


1.01 




29.54 


3.45 


N 


B2V,B1.5V 




W227 


140.3 


F 


1.58 


0.44 


38.61 


31.65 


3.45 


Y 


B2V,B1.5V 


Y 


W269 


16.7 


Q 


1.64 


0.62 




30.89 


3.45 


N 


B2V,B1.5V 




W228 


195.8 


F 


0.68 


4.70 




31.21 


3.32 


Y 


B2V 




W251 


50.45 


F 


0.35 


2.17 




30.44 


3.32 


Y 


B2V 




W607 


145.9 


n 

V 


0.43 


0.43 




31.04 


2.81 


N 


B3V 




W371 


225.6 


F 


0.44 


9.05 




31.23 


3.17 


Y 


B4V,B0.5V 




W336 


8.7 


Q 


0.12 


0.3 




29.40 


3.32 


? 


B5m 




W276 


93.76 


F 


0.25 


2.93 




30.71 


2.14 


N 


B6V 




W364 


7.6 


Q 


0.51 


2.25 
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*: F spectral fit, Q quantiles, C conversion from the count rate. 
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*** : ? no light curve. 
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Fig. 10. — log (Lx) vs. log (Ltoi) for the detected O (crosses) and B (circles) 
stars in NGC 661 1. The dashed line is the best fit line calculated for for stars 
with log(Ltoi) > 37 erg ■ s -1 , the solid line is the typical Lx ~ 10~ 7 ■ Ltoi 
relation. 

emission from O stars has been observed: a moderate hard 
component (kT = 2 - 3 keV) rotation ally modulated or with 
rapid variability (IStelzer et al.l 120051) . or a very hard ther- 
mal component due to shocks in the magnetically confined 
wind in the stellar equatorial plane in s tars with strong mag- 
netic dipole (Babel & Montmerle, 1997), wind-wind collision 
in close massive binaries (Pol lock et al.L 120051). which ca n be 
also source of soft x-rays emission dGagne etaTL I20TT1) . or 
a hard non-thermal component caused by inverse Compton 
scattering of UV stellar photons by relativisti c charged parti- 
cles in the stellar wind dChen & W hite. 199TJ). 

As listed in Table Q] all the O stars in NGC 661 1 have soft 
spectra (with one exception that will be discussed later), sug- 
gesting that in M16 the shocked wind emission is the only 
mechanism producing X-rays in the O stars. Fig. [I0]shows the 
log (Lx) vs. log (Lboi) relation for the OB stars in NGC 661 1. 
The best linear fit for stars with Lboi > 10 37 erg ■ s~ l (the 
sample of B stars with log (Lboi) < 37 erg ■ s~ l is strongly 
incomplete) gives a slope of 0.7 x 10~ 7 (dashed line), simi- 
lar to the L x ~ 10~ 7 • Lboi relation (represented by the solid 
line) typical o f the sh o cked wind emiss i on, as found first by 
lHarnden et all (119791) : iPallavicini et al.1 dl981l) and then ob- 
served in o ther massive star-f orming regions (i.e. in the Ca- 
rina region, INaze et ai.Ll20TTl) . 

A significant level of hard emission should be expected 
in principle from the 04V and 07.5V binary system W205, 
from the collision of the intense stellar wind emitted from the 
two components (expected to be ~ 2.0 x 10~ 6 M Q /year and 
~ 0.35 x 10~ 6 Mn Near: iSmithl [20061). but the sep aration be- 
tween the two stars is too large dSana et al.L 120091 reported a 
rotation period of hundreds of years) to expect a significant 
wind+wind X-ray emission. 

Only for W175 the spectral fit suggests the presence of a 
possible hard component. This system is composed of an 
08.5V and an 05V stars, with a third intermediate mass 
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Fig. 11. — Thermal plasma model for W175. The two temperature compo- 
nents are shown (dotted lines). 



(spect ral type A-F) component reported by Duche ne et al.l 
d2001l) . The separation bet ween the two mor e massive com- 
ponents is about 1200 AU (ISana et all 120091) . As expected, 
W175 is very bright in our Chandra observations, with 1258 
net counts and log (Lx) = 32.65 erg ■ s -1 . The spectral fit with 
abundances fixed at standard values is very poor, with xl = 
2.06 indicating a statistically unacceptable fit. An acceptable 
spectral fit (xl = 0.86) was instead obtained adopting a 2 tem- 
perature thermal model with abundances for Iron and Oxygen 
allowed to vary. The values of Nh, kT and Lx of W 175 ob- 
tained with this fit are reported in Table Q] and the abundances 
found are O = 2.3+0.7 Q and Fe = 0.7+0.2 Fe Q . In O stars a 
subsolar abundance of iron can be ex pected, as predicted, fo r 
example, by the model developed by IZhekov & P alla (2007), 
but in the case of W 175 the iron abundance is only marginally 
subsolar. This spectral fit suggests the presence of a hard com- 
ponent with a temperature kT = 4 ± 2 keV, even if the emis- 
sion measure of this hot plasma is more than 20 times smaller 
than that of the cold component, that can be produced in the 
wind-colliding zo ne such as ob s erved in other similar systems 
(i.e. HD 93 129 A. iGagne et all |20T1 . Fig. E] shows the 2T 
thermal plasma model fitted to W 175 spectrum. 

5.2. X-ray emission from B and intermediate-mass stars 

In the sample of early-B stars of NGC 6611 the ratio 
of bolometric flux emitted in X-ray ranges from -5.4 < 
log (Lx/ Lboi) ^ -8.5. This large scatter is related to the va- 
riety of mechanisms for X-ray production operating in this 
family of stars. Among these B stars, 10 have kT < IkeV 
and constant light curve, according to the KS test performed 
by AE. Their X-ray emission is then typical of shocked wind 
emission. On the other hand, 21 early B-stars have a plasma 
temperature kT > 1 keV, with six of them showing very hard 
spectra, with kT > 3 keV. For the B stars with plasma tem- 
perature kT > 1 keV the KS test cannot help to discern be- 
tween the possibility that the hard emission is due to an un- 
resolved low-mass companion or to the action of one of the 
mechanisms described above. However, among the 6 stars 
with kT > 3 keV, five of them are not spectroscopically classi- 
fied as close binary systems, and then for them the hard X-ray 
emission cannot be produced by colliding winds. Their light 



13 



curves produced by AE suggest some flaring activity that will 
be studied in detail in a forthcoming paper. This emission can 
come from a low-mass unresolved companion in a wide orbit 
observed during a flaring a ctivity, or from the B stars them- 
selves, since, as proposed by Feigel son et al.l (120021) . magnetic 
reconnection events can take place on the surface of young B 
stars. 

One peculiar case is W601. This B1.5V star has been sug- 
gested to be very young (0.016 Myr) and in a transition phase 
between the Class II and the MS phases, with a P Cygni- 
like H n line typ i cal of the presence of outflow/infall activity 
(Mart avan et all 120081) . This star has been studied in detail 
bv lAlecian et al.l d2008l) . who found an accretion rate equal to 
1 0~ 4 M -y r~ 1 , a rotational velocity of 1 90 + 1 Km -s~ l , and an 
intense mean magnetic field equal to 3 kG, classifying it as a 
Herbig Be star. With these characteristics, it is not surprising 
that we detected an intermediately hard component in W601 
spectrum (kT = 1.72 fee V). The KS probability that the light 
curve is constant is high (0.96). 

As expected, the fraction of detected sources in the A to 
late-B range of spectral classes is very low: 19% of late-B 
stars and 10% of the A-F stars, with log (Lx) < 30.4 erg ■ 
consistent with emission from a low-mass companion follow- 
ing the Lx vs. mass relation shown in Fig. [8] One evidence 
that the X-ray emission detected from A-late B stars is due 
to a lo w mass companion has been obtained bv lDe Rosa et alj 
d201 ll) . who found that multiplicity in a sample of X-ray de- 
tected intermediate-mass stars is about four times higher than 
in a control sample of stars. 

6. X-RAY PROPERTIES OF THE POPULATION OF 
SELECTED REGIONS OF Ml 6 

In this appendix we analyze the X-ray properties of the 
YSO populations of selected regions of M16 with different 
characteristics. Fig. [T2l shows these regions on an IRAC [8.0] 
image. In the following discussion, the me dian age of the 
members falling in each region are taken from lGuarcello et akl 
(2010b), which also includes a detailed description of the evo- 
lutionary status of the sources falling in some of these regions. 

6.1. The central cavity 

The central cavity contains the core of NGC 661 1, with al- 
most all the OB stars discussed in Sect. [5] and a rich YSO 
population (385 members detected in X-rays, among which 
21% have a disk) with a median age of lMyr. This region 
has the highest fraction of members with disks among all the 
regions analyzed in this section. The median Nh of the X- 
ray members is 0.52 x 10 22 cm~ 2 , corresponding to an ex- 
tin ction of Ay = 2.91 m , slightly larger than that extimated 
by iGuarcello et al.l (|2007j) from the optical color-magnitude 
diagrams (Ay = 2.7 m ); the median plasma temperature is 
kT = 1.95 keV, lower than the median value of the whole 
population. In the central cavity there is a significant popu- 
lation of candidate background sources (45 sources, the 10% 
of the total number of X-ray sources in this field). It is likely 
that most of them are really background contaminants, since 
we do not expect a large embedded population in this field 
where the molecular cloud has been cleared by the energetic 
radiation from the OB stars lying in this region. 

6.2. The pillars of creation 

The members of the pillars of creation are 33 (8 with 
disks), their age is about 2-3 Myr old in the bottom of the 




Fig. 12. — IRAC [8.0] image of M16 with marked all the regions whose 
population is analyzed in Sect. [6] 

pillar, in the south-east, and it gets younger in the direc- 
tion of the tip of the pillars. The median absorption of the 
X-ray sources falling in this region is typical of the cluster 
(Nh = 0.45 x 10 22 cto~ 2 ), suggesting that almost all these 
sources belong to the cluster but they are not physically as- 
sociated with the pillars. The median plasma temperature 
is kT = 2.03 keV, smaller than the median of the mem- 
bers of the whole cluster. Three X-ray sources are associ- 
ated with the tip of the largest pillar, called Pillar 1, where 
a protostar with a mass o f 4 - 5 M R and an intense infrare d 
excesses has been found (McCaughrean & Andersen, 2002). 
One (at a = 18 : 18 : 50.88 and 6 = -13 : 48 : 43.48) 
has been classified as a foreground source, and its Nh, de- 
rived from the quantile analysis, is compatible with this clas- 
sification (Nh = 2 x 10 19 cm~ 2 , corresponding to an Ay = 
0.01"'). The second source is a disk-less member located 
at a = 18 : 18 : 49.78 and 5 = -13 : 48 : 56.2. It 
has a spectrum fitted with a IT thermal plasma model, with 
kT = 1.71 feV and N H = 0.43 x 10 22 cm" 2 , which corre- 
sponds to the average extinction of the cluster suggesting that 
this source is not embedded in the pillar. The only likely X- 
ray counterpart of the protostar in the tip of the pillar is at 
a = 18 : 18 : 50.31 and 6 = -13 : 48 : 54.31, with 
103 net counts, an absorption Nh = 3.89 x 10 22 cm~ 2 , cor- 
responding to an Ay = 21.7™, and a plasma temperature of 
kT = 10.89 keV and log(Lx) = 31.69 erg ■ s^ 1 suggesting an 
embedded very active young star with a X-ray thermal spec- 
trum. The X -ray counterpart of this protostar has already been 
identified by iLinskv et alj (12007b . 

6.3. The Column V 
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In the other pillar of M16, Column V, fall 32 members (6 
with a disk) with an age of about 1 Myr. Their median absorp- 
tion is Nh - 0.95 x 10 22 cmr 2 , twice that of the central cavity, 
with the most obscured sources at the bottom of the pillar. The 
median value for the plas ma temperature is kT = 2.18 keV. 
Meab urn & Walshl (Q986) discovered a bright Herbig-Haro 
object in the tip of Column V. We detected one X-ray source 
at less than 1" away from it, at a = 18 : 19 : 04.65 and 
6 — -13 : 45 : 32.31. The quantile analysis revealed a 
plasma temperature of kT = 3.4 ke V, an intermediate ab- 
sorption (Nh = 0.3 x 10 22 cm~ 2 ) and an X-ray luminosity 
Lx = 10 erg ■ s -1 . Another candidate embedded proto- 
star of Column V detected in X-rays at a = 18 : 19 : 07.08 
and 6 - -13 : 45 : 22.54 can be associated with water masers 
identified bv lHealv et alj (120041) . This source is without a stel- 
lar counterparts in our catalog. It has been observed both in 
the c - field (32.7 net counts, log(Lx) = 30.6 erg ■ s _1 , and 
a spectrum fitted by a power-law) and the e - field (11.8 net 
counts). 

6.4. The Bright Rimmed Cloud SFO30 

Another region rich in cluster members is the Bright 
Rimmed Cloud SFO30, with 85 members (6 with a disk) de- 
tected in X-rays. In the optical and infrared diagrams, these 
sources are largely extincted (Ay > 5'"). The median ab- 
sorption (Nh = 1.16 x 10 22 cm~ 2 ) is similar to that of the 
sources falling in Column V, with a median plasma temper- 
ature of kT = 2.67 keV. The distribution of plasma tem- 
perature of these stars, however, has a very hot tail, with 
the 66% terzile at kT = 5.33 keV. Considering all the X- 
ray sources falling in this region, there is the largest frac- 
tion of candidate members without a disk and the lowest of 
can didate background sources. Since, from infrar ed stud- 
ies dlndebetou w et all 120071: iGuarcello et all 12009ft and de- 
tection of water masers (Hea lv et all, [2004), a population of 
embedded protostars is expected to fall in this region, it is 
possible that the extinction here is high enough to efficiently 
absorb all the background sources or that these protostars are 
not X-ray active yet (the population of this region is on aver- 
age younger than 1 Myr). 

6.5. The north-east embedded cluster 

The most extinguished region is the N E molecular cloud, 
where a young embedded cluster lies dlndebetouw et all 
2007; Guarcellp_etali|2009). This region has the highest frac- 
tion of X-ray sources without stellar counterparts (82 sources, 
52% of all X-ray sources falling in this region). Among the 70 
candidate YSOs in this region, 12 have a disk and are classi- 
fied as embedded YSOs from their IRAC colors (among them 
also the Class I sources detected in X-ray discussed in Sect. 
0). It is not surprising that this is the region with the high- 
est median absorption, with Nh = 1.92 x 10 22 cm~ 2 , which 
is four times larger than the value in the central cavity and 
corresponds to an Ay = 10.7 m . 

6.6. The other regions 

The remaining regions are the poorest in candidate cluster 
members. In the East region the 60 members detected in X- 
rays (12 with a disk) are older than 1 Myr, with a median ab- 
sorption similar to the central cavity (Nh = 0.65 x 10 22 cmr 1 ) 
and a median plasma temperature of kT - 2.5 keV. The 
members detected in X-ray falling in the West region (25 
sources, none with a disk) have similar ages, an absorption 



lower than the members in the central cavity (median value 
Nh - 0.33 x 10 22 cm~ 2 ) and a median plasma temperature 
of kT = 2.45 keV. The South-East region contains the oldest 
X-ray detected cluster members (12 sources, one with a disk, 
with a median age of about 3 Myr), with the lowest median 
plasma temperature (kT - 1 .72 keV), and a median absorp- 
tion of the candidate members detected in X-rays smaller than 
that of the central cavity (N H = 0.36 x 10 22 cm" 1 ). 

7. SUMMARY AND CONCLUSIONS 

In this paper we analyze three 80 ksec Chandra/ACIS-I ob- 
servations of the Eagle Nebula: one archival centered on the 
young open cluster NGC 6611 in the center of the nebula, 
one centered on the evaporating pillars named "Column V", 
and one in the north-east region centered on an embedded 
cluster of very young stellar objects. The final catalog of 
the X-ray sources amounts to 1755 entries, classified from 
their photometric properties as "disk-bearing cluster mem- 
bers" (219 sources), "disk-less cluster members''' (964), "fore- 
ground sources" (76), and "background sources" (504). 

Spectral properties and X-ray luminosities have been de- 
rived from spectral fitting and analysis of the photon energy 
quantiles. Considering all the cluster members detected in 
X-ray, we found a median Nh = 0.640 x 10 22 cm~ 2 , corre- 
sponding to a median extinction of Ay = 3.6™, and a me- 
dian kT = 2.13 keV. The median Nh is three times larger 
in the North-East field with respect to the central field where 
NGC 661 1 has cleared most of the parental cloud. 

We compared the X-ray properties of M16 members and 
those in the Orion Nebula Cluster, which is the best charac- 
terized sample of X-ray sources in a young cluster. The X-ray 
Luminosity Function has similar slope (a power law with T = 
-0.85 +0.09) in a range of luminosities log (Lx) > 30erg-s~ l , 
which is the completeness limit of our observations. This al- 
lowed us to estimate that the total population of M16 should 
be about 2700 members. Also the slope of the Lx vs. mass 
function is similar to that in Orion for masses smaller than 
0.8 Mq, while for larger masses the distribution is flat around 
log(Lx) - 30.72 erg ■ s~ x , with a large spread. These analo- 
gies support the universality of the X-ray luminosity functions 
and the L x vs. mass relation for clusters younger than 5 Myrs. 

Our study supports the evidence that the X-ray activity in 
disk-bearing stars is less intense than that in disk-less stars 
of similar mass. This result has been obtained also in other 
young clusters, while in some other cases no dependence with 
the presence of disk in the X-ray activity has been observed. 
The present study of M16 supports then the idea that the pres- 
ence of a circumstellar disk, likely due to the accretion pro- 
cess, affects the coronal activity in the hosting star. 

We took advantage of the large population of massive stars 
detected in X-rays (11 O stars, a 85% detection rate, and 31 
B stars, 39%) to study the mechanism for the production of 
X-rays in such massive stars. All the O stars but one have 
soft spectra without a hard component and have constant light 
curves, suggesting that the X-ray photons are produced by 
shocks propagating in the stellar wind. A hard component 
has been detected only in one case (a triple system with two 
O stars and a A-F star), but it is only marginally significant, 
contributing less than the 5% of the total spectrum. 

The sample of detected B stars is more heterogeneous, with 
10 sources with soft spectra typical of wind emission, and 21 
with kT > 1 keV. For these latter sources it not possible to 
discern the possibility that some mechanism for the produc- 
tion of hard X-ray photons is at work or the emission is due to 
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unresolved companions. However, the binarity of the massive 
stars of Ml 6 has been studied by several authors and most 
of these stars do not show evidence for such companions. In 
particular, among the six B stars with kT > 3 keV, only one 
is a binary system, but with 2 B stars, and all of them seem 
to have flares in their light curves. This result suggests that 
that in these young massive stars some magnetic reconnection 
events can occur on the stellar surface. 

Finally, the X-ray properties of the sources in relevant re- 
gions of Ml 6, such as the central cavity and the Pillars of Cre- 
ation, have been discussed to analyze the different properties 



of the M16 population in different locations of the nebula. 
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APPENDIX 
SOURCE CATALOG 

The information about the X-ray sources of M16 analyzed in this work are available in an on line catalog. The catalog is not 
shown here because of the large number of columns. In the following we describe the data available in the catalog. 

• RA, DEC: source position in J2000. 

• <Tpos'- the uncertainty on source position found by AE. 

• ft the source off-axis angle. 

• PBS : the probability associated with the null-hypothesis that there is not a real source associated with the position. Only 
21 sources have the PBS value larger than the threshold of 0.1 typically used with AE. 

• Sgnf: photometric significance computed as the ratio between the net counts in 0.5 - %keV on the upper error of net 
counts. 

• Status: class of the source, based on its optical/infrared properties (disk-bearing members, disk-less members, foreground 
source, background source). 

• Cnts: the total number of counts associated with the source in the whole energy band (0.5 - 8 keV). 

• Bkgjonts: the observed background counts associated with the source in the whole energy band (0.5 - 8 keV). 

• Net-cnts: the net counts associated with the source in the whole energy band (0.5 - 8 keV). 

• PS F_frac: the fraction of the PSF at 1 .497 ke V corresponding to the extraction region. 

• medE: the median photon energy in keV. 

• Nh'- the value of the hydrogen column density in the direction of the source, in units of 10 22 crrf 2 . This value is estimated 
either from the spectral fitting or quantile analysis as described in Sect. [3] with errors. 

• kT: temperature of the emitting plasma, in keV. This value is estimated either from the spectral fitting or quantile analysis 
as described in Sect. [3] with errors. 

• kT2'- second temperature of the emitting plasma, in keV, for the sources whose spectra are well-fitted by a thermal 2- 
temperature plasma model, with errors. 

• Pq: index of the power law for the sources whose spectrum is best fitted by a power-law model. 

• Lx'- the adopted X-ray luminosity for the source. See sect. |4]for details on how Lx is determined, with errors. 

• Model: the model used to fit the source spectrum (thermal, 2T thermal, power law). See Sect. 12.31 

• Pks'- probability based on the K-S statistic with the null hypothesis that the source is not variable (in the total band). The 
classification suggested by the AE manual is: no evidence for variability (0.05 < Pks)\ possibly variable (0.005 < Pks < 
0.05); definitely variable (P KS < 0.005). 

All the information, but the position, PBS and Status, are provided for all the three observations, sorted in alphabetical order 
(c - field, e - field, and ne - field), for a total of 63 columns. 
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